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Abstract 

Light sterile neutrinos with masses at the sub-eV or eV scale are hinted by cur- 
rent experimental and cosmological data. Assuming the Major ana nature of these 
hypothetical particles, we discuss their effects in the neutrinoless double beta de- 
cay by exploring the implications of a vanishing effective Majorana neutrino mass 
(m)ee- Allowed ranges of neutrino masses, mixing angles and Majorana CP-violating 
phases are illustrated in some instructive cases for both normal and inverted mass 
hierarchies of three active neutrinos. 
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1 Introduction 



Despite the success of the standard three-neutrino oscillations in explaining the results 
of solar, atmospheric, reactor and accelerator neutrino experiments PQ by two distinct 
mass-squared differences, a number of anomalies from the short base-line antineutrino 
experiments (e.g., LSND [2] and MiniBooNE [3]) and the latest reactor antineutrino 
fluxes [1] indicate the existence of oscillations with much shorter baselines. This would 
imply the existence of extra mass-squared differences and therefore the mixing of active 
neutrinos with extra sterile neutrino states [51 [6]. Furthermore, the analysis [7] of cosmic 
microwave background and large scale structure data is hinting an existence of additional 
radiation in the Universe, with sterile neutrinos being one of the plausible candidates. 
Finally, one extra sterile neutrino is also allowed and even favored from a recent analysis 
of the Big Bang Nucleosynthesis bound [8]. Therefore, we should be open-minded on their 
existence and it might be instructive to study the light sterile neutrino hypothesis with 
other oscillation channels [9] or in non-oscillation processes including the beta decay [10] 
and neutrinoless double beta decay [TT1 [12] . 

The most promising process to probe the Majorana nature of massive neutrinos and 
acquire the information of Majorana CP-violating phases is the neutrinoless double beta 
decay (Qv(3f3) [121 US] of some even-even nuclei: 



which takes place via the exchange of Majorana neutrinos but would be prohibited if 
neutrinos are the Dirac particles. The decay rate is proportional to an effective Majorana 
neutrino mass term (m) ee and to the associated nuclear matrix element. The latter can 
in principle be calculated with some uncertainties [32] in nuclear physics. The effective 
neutrino mass (m) ee is related to the neutrino masses, mixing angles and Majorana CP 
phases as follows: 



where rrii is the mass of the i neutrino mass state and the sum of i includes all light 
neutrino mass states. V e i stands for an element in the first row of the lepton flavor mixing 
matrix, and it contains a Majorana CP phase in general. 

Great efforts have been made to search for signals of the 0^/3/3 process [T2], [13] recently, 
however, only negative results have been obtained and the upper bound on the effective 
neutrino mass is achieved accordingly (i.e., (m) ee < 0.3 — 1.0 eV [33]). Although a positive 
Ov/3/3 signal would verify the Majorana nature of massive neutrinos, negative results do 
not necessarily mean that neutrinos are the Dirac particles. For Majorana neutrinos, a 
vanishing effective neutrino mass (m) ee due to cancelation among different mass eigen- 
states could also lead to the failure of Qvf3/3 observations. This possibility in the standard 
scenario of three- neutrino mixing has been explored in some previous works [TU [T7] . If 
light (sub-eV) sterile neutrinos exist and are mixed with their active counterparts, they 
will also contribute to the 0^/3/3 decay rate and modify the prediction on the effective 
neutrino mass (m) ee . Therefore, the parameter space of a vanishing {m) ee will get altered 
drastically in contrast to the standard scenario. This point constitutes the main concern 
of the present work. We shall derive the correlative relations among neutrino masses, 






(2) 
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mixing angles and Majorana CP phases and constrain the corresponding parameters by 
using current experimental data of active and sterile neutrino oscillations. 

The main part of this work is organized as follows. In section 2, we shall describe the 
theoretical framework and derive the analytical formulas for a vanishing (m) ee . Section 3 
is devoted to a numerical analysis. The analytical relations will be confronted with current 
global oscillation data and numerical examples will be illustrated in several specific cases. 
Finally, we conclude in section 4. 



2 Analytical Calculations 

Assuming that neutrinos are the Majorana particles and there exist N s species of light 
sterile neutrinos, the neutrino flavor eigenstates (u e , v T ) are connected with the 3 + N s 
mass eigenstates (z/«, i — 1, ... ,3 + N s ) by a 3 x (3 + N s ) flavor mixing matrix V in the 
charged- current interactions, which contains 3 + 3N S mixing angles, 1 + 2N S Dirac CP 
phases and 2 + N s Majorana CP phases. In the 0^/3/3 process, as revealed in Eq. (2), 
only the mixing matrix elements in the first row of V are relevant in our calculations. 
Without loss of generality, we may redefine the phases of three charged lepton fields in 
an appropriate way such that the phase of the mixing matrix element V ei is only of the 
Majorana type. So the expression for the effective neutrino mass (m) ee can be rewritten 
as 



(m) e 



3+N s 

m 1 \V el \ 2 e 2ip ^ + m 2 \V e2 \ 2 e 2ip * + m 3 \V e3 \ 2 + £ mj \V ej \ 2 e 2ip 



j=4 



(3) 



where p i (i — 1, . . . , 3 + N s ) are the Majorana phases and a vanishing p 3 has been chosen. 
Although some schemes with more than one sterile neutrino may have their distinct 
properties such as CP violating effects in active-sterile neutrino oscillations [3 E] , we can 
combine all the sterile neutrino contributions in the 0uj3/3 process by assuming an effective 
(3 + 1) scheme in which the additional neutrino mass m , mixing matrix element \V e0 \ 
and Majorana CP phase p are defined by 

3+N s 

m \V e0 \ 2 e 2ip o = Y, mj y,j-<- : "- ■ (4) 

j=4 

Taking account of Eq. (3) and (4), we find that a vanishing effective neutrino mass 
(i.e.,(m) ee = 0) requires 

m \ V e o\ 2 sin2p + mj V e i\ 2 sin 2p l + m 2 \ V^] 2 sin 2p 2 = 0, 
^o\V e o\ 2 cos2p + m^Vei] 2 cos2p 1 + m 2 \V e2 \ 2 cos2p 2 + m 3 \V e3 \ 2 = 0. (5) 

Comparing these two conditions with current experimental data on the mixing matrix 
elements (|V el |, \V e2 \, \V e3 \ and |V^ |) and the mass-squared differences (Am 21 = m 2 — m 2 , 
Amgj = m 2 — mf and Am§j = tuq — m 2 ), we might be allowed to determine or constrain 
the absolute neutrino mass scale, neutrino mass hierarchies and Majorana CP phases. 
Before doing the most general analysis, however, let us discuss several specific cases to 
reveal the main consequences of Eq. (5). 
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(a) The case of m = or \V £ q\ = (or both) means that contributions to (m) ee are 
canceled out among different sterile species of the mass eigenstates, which is practically 
equivalent to the standard scenario of three-neutrino mixing. Taking the (3 + 2) scheme 
with two additional sterile neutrinos as an example, the cancelation between v A and u 5 
indicates 



m 4 \V e5 \ 2 2n + l 

m 5 =~WJ^ l*-*l = — (6) 



with n being an arbitrary non-negative integer. In this vanishing (m) ee only 

requires the parameter correlations among three active neutrinos, 



mi 
m 2 





v e2 \ 


2 sin 2p 2 ni2 


v e3 \ 


2 sin2pi 




v el \ 


2 sin 2pi ' m 3 


v e2 \ 


2 sin (2p 2 - 2pi) 



(7) 



which have already been studied in some earlier literatures [T6| ITT] . Therefore, we shall 
only discuss the process with non-vanishing sterile neutrino contributions in the 

following parts of this work. 

(b) The case of CP invariance requires that all the Majorana CP phases should take 
some specific values (i.e., p i = niir/2 with rij being arbitrary integers). So the two 
conditions in Eq. (5) will be simplified to a single one, 



-l)'°m |Ko| 2 + (-l^ilKil 2 + (-l) l2 m 2 \V e2 \ 2 + m 3 \V e3 \ 2 = 



with l k = or 1 for (k = 0, 1, 2), depending on n k being even or odd integers respectively. 
Thus the mass spectrum of active and sterile neutrinos might be fully determined for each 
group of (Iq, li, l 2 ) with the help of current knowledge on neutrino mass-squared differences 
and mixing matrix elements. A simple calculation shows that (0, 1, 0) and (1, 0, 1) are 
permitted for both mass hierarchies of active neutrinos, (1, 0, 0) and (0, 1, 1) are allowed 
only for the case of Am^ > 0, but all the other possibilities are ruled out by current 
global neutrino data. 

(c) The case of a neutrino mass eigenstate being massless depends on the sign of Am^. 
m 1 = or m 3 = might hold for normal or inverted mass hierarchies respectively. If 
Amjj > and m 1 = hold, then Eq. (5) requires 





V e0 \ 


2 sin2p m 3 _ 


v e0 \ 


2 sin (2p - 2p 2 ) 




v e2 \ 


2 sin 2p 2 ' niQ 


v e3 \ 


2 sin 2p 2 



m 2 
niQ 

on the other hand, in the case of Am| x < and m 3 = 0, we are led to 

m 2 
m 



(9) 



\V e0 \ 


2 sin (2p 


-2pi) 


nil 


\VeO 


| 2 sin (2p 


-2p 2 ) 


\v e2 \ 


2 sin (2p 2 


- 2pi) ' 
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\V e l 


2 sin (2p 2 


-2pi) 



(10) 



Since the masses of active neutrinos are totally calculable when we assume m 1 = 
or m 3 = 0, part of the parameter space of Majorana CP phases must be excluded. 
Furthermore, the allowed region of (m , |V e o|) might be constrained either by using Eq. 
(9) or Eq. (10) from the active neutrino data or by using Eq. (4) from the sterile 
neutrino data. Therefore, a comparison between the two sectors may test the conditions 
of (m) ee = and give complete constraints on the (m , |14o|) parameter space. 
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In addition to the special cases considered above, we can derive the most general 
results about (m , \ V e o\) and p from the conditions in Eq. (5). Indeed, 

m o\ V eo\ 4 = mi\V e i\ 4 + ml\V e2 \ 4 + ml\V e3 \ 4 + 2ra 1 m 2 | V el \ 2 \ V e2 \ 2 cos (2p 1 - 2p 2 ) 

+2m 1 m 3 |Ki| 2 |K3| 2 cos2pi + 2m 2 m 3 | Ve 2 | 2 | V e 3| 2 cos2p 2 , (11) 

as well as 



m l\ 




2 sin 2pi 


+ m 2 


V e2 | 


2 sin 2p 2 
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2 cos2pi + m 2 \ 


v e2 | 
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cos2p 2 + m 3 \ 
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(12) 



Note from Eq. (11) and Eq. (12) that the parameters of sterile neutrinos (i.e., ?n |Veo| 2 
and po) may be calculated from the active neutrino data, depending on the absolute 
neutrino mass scale and neutrino mass hierarchies. On the other hand, we may directly 
constrain m^Veo! 2 an d Po from the sterile neutrino data by using Eq. (4). A combined 
analysis of both sectors will give the complete parameter space allowed by the global 
neutrino data. 



3 Numerical Analysis 

To be explicit, let us carry out a numerical analysis of the implications of (m) ee = 0. 
For illustration, we shall take the (3 + 2) neutrino mixing scheme with two sub-eV 
sterile neutrinos as an example, which is more or less favored by the latest experimental 
[5, 6] and cosmo logical data [7J. We may also assume that these sterile neutrinos under 
consideration do not significantly affect the values of two mass-squared differences and 
three mixing angles of active neutrinos extracted from current active neutrino data [IS] 119). 
In this assumption we shall use the latest results on active neutrino parameters in Ref. 
[j~8] and sterile neutrino data in Ref. [5], namely, Am 2 x « 7.58 x 1(T 5 eV 2 , Am 2 x « 
2.39 (-2.31) x 10" 3 eV 2 , Am^ « 0.47 eV 2 and Am 2 51 rj 0.87 eV 2 together with |V el | « 
0.824, \V e2 \ « 0.547, |y e3 | « 0.145, \V e4 \ « 0.128 and |V; 5 | « 0.138 as typical inputs 
of the best-fit values. The sign of Am 2 x is unknown and the absolute neutrino mass 
scale, which can be characterized by the smallest neutrino mass {m 1 for Am 2 x > or 
m 3 for Amjj < 0), is expected to be at most of O(0.1) eV or 0(1) eV as constrained 
by current experimental and cosmological data [U[20]. On the other hand, the Majorana 
CP phases are totally unconstrained and will vary in the full parameter space. Finally, 
we shall take into account the uncertainties of active neutrino parameters at the la or 3cr 
confident level in our numerical analysis but neglect the uncertainties of sterile neutrino 
parameters which are not explicitly presented in Ref. [5] . We admit that our calculations 
are quite preliminary and mainly serve for illustration, but it should be emphasized that 
they indeed reveal the main consequences of (m) ee = 0. 

Now we are able to do the numerical analysis on possible implications of (m) ee = 0. 
We shall employ the analytical results of different cases discussed in section 2 as our 
guidelines and plot the corresponding regions allowed by the global oscillation data. 

In the cases of the smallest neutrino mass being zero, it is straightforward to present 
the correlative relations among different Majorana CP phases and between (m Q , |V e o|). 
First we illustrate the case of m 1 = in Fig. 1. The parameter space of two Majorana 

1 An analysis of the latest global data has also been done by a few independent groups, such as Ref. 
[T§] for the standard three-neutrino oscillations and Ref. [5] for the active-sterile neutrino oscillations . 
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phases (p , p 2 ) is shown in Fig. 1(a), where different (black and grey) colors stand for 
lcr and 3a ranges of the active neutrino parameters respectively. One should note that 
p 1 is irrelevant in this case because the massless mass eigenstate v\ does not contribute 
to the 0^/9/3 process. In Fig 1(a), we can learn that there are two isolated regions to 
guarantee (m) ee = 0, and they are related through the transformation of (p — > n — p , 
p 2 — > 7i — p 2 )- This point can be understood by considering the fact that the formulas in 
Eq. (9) are invariant under the same transformation. Furthermore, the allowed regions 
can be generalized to the larger (p + n 7r, p 2 + n2^) parameter space, where n and n 2 
are arbitrary integers. The allowed regions of (m , |Ko|) can be found in Fig. 1(b), where 
the (black and grey) shaded regions stand for constraints from the active neutrino data 
and the region with (red) sparse lines is required by the sterile neutrino parameters. The 
(black) horizontal line stands for the 2a upper bound on the sum of the neutrino masses 
from cosmological probes [2TJ and only the region below this line is allowed. The profiles 
of these allowed regions have the shape of the inverse-rate curve in the single logarithmic 
scale, which can be easily understood by the analytical expressions in Eq. (4) and Eq. 
(5). One should also note that there is no overlap between the allowed regions of active 
and sterile neutrino data even within the 3a ranges, indicating that the case of m 1 = 
and Awijj > is disfavored by current global oscillation data if (m) ee = holds. This 
conclusion demonstrates the fact that the minimum of m |Ko| 2 from the sterile neutrino 
data, | a/ Am^ | V^4 1 2 — a/ Am^l Ks| 2 | is always larger than the maximum of m |V^o| 2 from 
the active neutrino data, | A/Am|j V e2 \ 2 + \^Km^\V e3 \ 2 \. 

A similar analysis for Am^ < and m 3 = is given in Fig. 2. The allowed regions 
of Majorana CP phases and (m , | V^o I ) are illustrated in the upper and lower panels 
respectively. Only the differences of Majorana phases can be constrained because it 
is the third mass eigenstate 1/3 that decouples from the effective mass (m) ee due to its 
masslessness. In Fig. 2(a), we take (p — Pi) and (p 2 — Pi) as free parameters and find that 
they are strongly correlated and constrained. The area of the allowed region here is much 
smaller than that in Fig. 1(a) simply for the reason that the mixing matrix elements of 
active neutrinos in Eq. (10) (i.e., \V e \\, \V e 2\) are much more precise than that appearing 
only in Eq. (9) (i.e., \Ve^\). Now let us take a look at Fig. 2(b) for the allowed regions 
of (m , |V^o|)- The (black) horizontal line stands for the 2a upper bound on the sum of 
the neutrino masses from cosmological probes |21j. Its most distinct difference from Fig. 
1(b) is the overlaps between active and sterile neutrino constraints appearing at both la 
and 3a levels. Therefore, m 3 = and (m) ee = can be simultaneously satisfied in the 
presence of sub-eV sterile neutrinos. In the meantime, it is instructive to compare the 
present scenario with the standard scenario of three-neutrino mixing, where a vanishing 
(m) ee favors Am^ > over Am| x < 0. 

So far we have only discussed the consequences of a vanishing (m) ee in some specific 
cases. To illustrate the dependence on the absolute neutrino mass scale, we calculate 
the allowed ranges of m^V^I 2 versus the smallest neutrino mass (m 1 or m 3 ) for both 
neutrino mass hierarchies in Fig. 3. We give the constraints from the active neutrino 
data with (black and grey) shaded regions and constraints from the sterile neutrino data 
with regions of (red) sparse lines. The overlapping regions of both sectors are the complete 
parameter space of {m) ee = and the global oscillation data. The case of Am 31 > is 
illustrated in Fig. 3(a), which confirms the previous conclusion that the specific case 
of TTij — is disfavored by current global oscillation data. As uii is getting larger, the 
active neutrino spectrum changes from the hierarchical type to the quasi- degenerate one 
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and therefore the relations between m |V^o| 2 and m 1 in Eq. (11) approximate to linear 
functions. Numerically, we may find the allowed ranges explicitly as 0.005 eV < m 1 < 
0.116 eV and 0.006 eV < «\)l^eo| 2 ^ 0.029 eV. In comparison, a similar calculation for 
Amjj < can be found in Fig. 3(b) and the ranges of m 3 and m |V^o| 2 are < m 3 < 
0.10 eV and 0.012 eV < m |V e o| 2 < 0.029 eV, respectively. Finally, let us comment on the 
Majorana CP phase p in Eq. (12). Although it can be calculated from the active neutrino 
data, there is no correlation between p and the absolute mass scale (characterized by m 1 
or m 3 ). Our numerical calculation shows that p may run over the full parameter space 
at each point of the smallest neutrino mass (m 1 or m 3 ). Only when the neutrino mass 
scale is fixed, as the specific cases discussed above, there can be non-trivial correlations 
among different Majorana CP phases. Within the full parameter ranges, it is instructive 
to estimate the minimal value of the effective Majorana mass (i.e., (rn)™ n ) where {m) ee 
does not vanish. We calculate (jn)™ n versus the smallest neutrino mass (m 1 or m 3 ) for 
both Am| x > and Am| x < cases in Fig. 4. The black solid and grey dashed lines stand 
for la and 3a ranges of the active neutrino data respectively. When Am^ > holds, 
there are two separated ranges of mi for a non-vanishing (m)™ m . A hierarchical neutrino 
mass spectrum with mi < 10~ 3 eV predicts a minimal value of (m) ee as (3 ~ 4) x 10 -3 eV 
and the magnitude of (rn)™ n can reach the (9(0.1) eV level for a nearly degenerate 
spectrum of the active neutrinos. In contrast, a non- vanishing (m)™ n can only achieved 
for m 3 > 0.1 eV in the case of Am 31 < and (m)™ m grows as m 3 becomes larger. 

4 Concluding Remarks 

Motivated by non-observation of the 0v(3(3 process and possible hints of sub-eV sterile 
neutrinos from the latest experimental and cosmological data, we have discussed the 
implications of a vanishing effective Majorana mass (i.e., (m) ee = 0) in the presence of light 
sterile neutrinos. Current neutrino oscillation data from both active and sterile neutrinos 
do allow (m) ee = to hold if neutrinos are the Majorana particles and their masses, 
mixing angles and Majorana CP phases possess some specific correlations. The allowed 
parameter space can be drastically altered in contrast to that of the standard three- 
neutrino mixing scenario. For the standard scenario, a vanishing (m) ee favors Am 31 > 
over Am 31 < 0. However, it is not really the case after we include the light sterile neutrino 
contributions. In the specific case where the smallest neutrino mass vanishes, a vanishing 
(m) ee is consistent with the case of Am 31 < rather than Am 31 > 0. Furthermore, for 
the complete parameter ranges, both cases are allowed and we may constrain the smallest 
neutrino mass as 0.005 eV < m 1 < 0.116 eV and < m 3 < 0.10 eV respectively. 

Going beyond the Standard Model, there are different mechanisms [22] that may be 
able to mediate the process, including light (active and sterile) Majorana neutrinos, 
heavy sterile Majorana neutrinos [231 121] and other particles [22] that can induce the 
lepton-number-violating (LNV) interactions. We stress that we have simply neglected all 
the contributions other than light Majorana neutrinos and our conclusions will be invalid 
once other mechanisms can give a comparable decay rate as either light active neutrinos 
or light sterile neutrinos. In addition, when more accurate data on neutrino mass and 
mixing parameters and more sensitive observations or constraints on the 0z^/3/3 process are 
achieved in the future [26l [27] , it might be possible to confirm or rule out the parameter 
space we have obtained and even probe the existence of light sterile neutrino states. 
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Figure 1: The regions of (p , p 2 ) (upper panel) and (m , |Ko|) (lower panel) allowed 
by (m) ee = and current neutrino oscillation data with Am^ > and m 1 = 0. The 
black and grey scattered regions (upper panel) and shaded regions (lower panel) stand for 
la and 3a ranges of the active neutrino data respectively. The region with (red) sparse 
lines (lower panel) is constrained from the sterile neutrino data. The (black) horizontal 
line stands for the 2a upper bound on the sum of the neutrino masses from cosmological 
probes. 
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Figure 2: The regions of (p — p x , p 2 — Pi) (upper panel) and (m , |Ko|) (lower panel) 
allowed by (m) ee = and current neutrino oscillation data with Am^ < and m 3 = 0. 
The black and grey scattered regions (upper panel) and shaded regions (lower panel) 
stand for la and 3cr ranges of the active neutrino data respectively. The region with 
(red) sparse lines (lower panel) is constrained from the sterile neutrino data. The (black) 
horizontal line stands for the 2a upper bound on the sum of the neutrino masses from 
cosmological probes. 
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Figure 3: The regions of ^ l^eo| 2 versus the smallest neutrino mass (m 1 or m 3 ) with 
Amjj > (upper panel) or Am^ < (lower panel). The black and grey shaded regions 
stand for la and 3a ranges of the active neutrino data respectively. The region with (red) 
sparse lines is constrained from the sterile neutrino data. 



13 




10" 5 10~ 4 10" 3 10~ 2 10" 1 

mi (eV) 




rri3 (eV) 



Figure 4: The minimal values of the effective Majorana mass {{m)™ in ) versus the the 
smallest neutrino mass (m 1 or m 3 ) in the (3 + 2) scheme with Am^ > (upper panel) 
or Arajj < (lower panel). The black solid and grey dashed line stand for la and 3a 
ranges of the active neutrino data respectively. 
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